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ABSTRACT 

We report the discovery of a significant population of RR Lyrae variables in the dwarf spheroidal galaxy Leo I. 
Based on 40 V and 22 B images of the galaxy taken using the ESO Wide Field Imager we have identified so far 
74 candidate RR Lyrae's in two CCD's hosting the main body of the galaxy. Full coverage of the light variations 
and pulsation periods have been obtained for 54 of them, 47 of which are Bailey ab-type RR Lyrae's (RRab's) 
and 7 are c-type (RRc's). The period distribution of the presently confirmed sample of RRab's peaks at P=0?60, 
with a minimum period of 0?54. The pulsational properties indicate for Leo I an intermediate Oosterhoff type, 
similar to other dwarf galaxies in the Local Group and the LMC. However, the rather long minimum period 
of the ab-type variables, and the significant number of RRab's with long period and large amplitude, suggest 
that the bulk of the old population in Leo I is more like the Oosterhoff type II globular clusters. The most 
straightforward interpretation is that a range in metallicity is present among the RR Lyrae's of Leo I, with a 
significant population of very metal-poor stars. Alternatively, these OoII variables could be more evolved. The 
average apparent magnitude of the RR Lyrae's across the full cycle is (V(RR)) = 22.60 ±0.12 mag, yielding a 
distance modulus (m-M)v.o = 22.04 ± 0.14 mag for Leo I on the "long" distance scale. 

Subject headings: Galaxies: individual (Leo I) — galaxies: dwarf — galaxies: stellar content — Local Group — 
stars: horizontal-branch — stars: variables: other 



1. INTRODUCTION 

RR Lyrae have long been recognized to be excellent tracers 
of old stellar populations, as well as good distance indicators 
for Population II systems. RR Lyrae variables and/or extended 
horizontal branches (HB's) are now known to exist in all Local 
Group dwarf spheroidal (dSph) galaxies, despite the different 
and complex star formation histories (Mateo 1998; Da Costa 
1998). The Leo I dSph was believed to represent the only ex- 
ception to this general trend, with its predominantly "young" 
population, and a wealth of intermediate age He-burning stars 
giving rise to a conspicuous red clump (Lee et al. 1993; Ca- 
puto et al. 1998; Gallart et al. 1999a, 1999b). Recently, color- 
magnitude diagrams (CMD's) extended to the outer regions have 
revealed a horizontal branch well populated from blue to red 
(Held et al. 2000). Given the presence of an extended HB and 
the low metallicity of the system ([Fe/H]~ —2, Lee et al. 1993) 



one would naturally expect to find RR Lyrae variables in Leo I. 
Indeed, a few stars near the limit of photographic photometry 
were mentioned by Hodge & Wright (1978) as probable RR 
Lyrae's caught at maximum light. Unpublished work of Keane 
et al. (1993) also reports the presence of these variables. 

To confirm these early suggestions, a search for RR Lyrae 
and other short period variables in Leo I was undertaken using 
the CCD mosaic imager at the 2.2m ESO-MPI telescope. In 
this letter we report on the successful detection of a conspicu- 
ous number (> 70) of RR Lyrae variables in this galaxy. The 
identification of RR Lyrae's in Leo I and the determination of 
their average luminosity and pulsational properties provide us 
with important tools for studying the metallicity and age of the 
oldest stellar population and the early star formation history of 
Leo I, by allowing at the same time a new independent mea- 
surement of the distance to this galaxy. 



1 Based on data collected at E.S.O. La Silla, Chile, Prop. No. 65.N-0530 
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2. OBSERVATIONS AND REDUCTION 

Images of Leo I (a20oo=10 08 27, (52ooo=12 18 30) through 
the B, V, I filters were obtained on the nights of April 20-25, 
2000 with the Wide Field Imager (WFI) at the 2.2 m ESO-MPI 
telescope at La Silla, Chile. The WFI camera consists of 8 
"ccd 44"-type EEV CCD's with a total field of view 34' x 33'. 
This coverage allowed us to observe the body of the galaxy and 
its outer regions in just one pointing. The main body of the 
galaxy was centered on the lower four CCD's of the mosaic to 
reduce light contamination from Regulus, which is about 20 ar- 
cmin south. The full data-set comprises 22 B, 40 V, and 4 / 
15-min exposures totaling <~ 20 hours of observing time spread 
over 6 consecutive half nights. Weather conditions were gen- 
erally photometric, while the seeing varied considerably during 
the observations. Standard stars from Landolt (1992) were ob- 
served during all photometric nights of the run in all CCD's. 

Basic reduction (flat fielding, registration, co-adding) of the 
CCD mosaic data were performed using the IRAF 3 package 
MSCRED (Valdes 1998) and the pipeline script package WF- 
PRED developed at the Padua Observatory. Reductions were 
done, and data calibrated, for the individual CCD's, yielding 
calibration uncertainties of 0.03 mag in B and 0.04 mag in V 
(night-to-night rms scatter of the zero points). For variable 
stars, the calibration were applied only to average magnitudes 
and colors given the relatively large color terms of the WFI 
camera (about 0.28 and -0.075 in b and v, respectively). The 
ALLFRAME program (Stetson 1994) was used to measure all 
images and combine the individual photometric measurements 
into a master catalog for each CCD. The combined photometry 
reaches V ~ 25.5 (about 3 mag fainter than the HB) with errors 
oy = <tb = 0.2 mag (standard error of the mean). The limiting 
magnitude for the individual frames varies between V = 23.0 
and 24.5 depending on the seeing and photometric conditions. 
Measurement errors vary accordingly between oy = 0. 12 - 0.35 
at V <~ 23.5. The measurement errors for stars at the HB level, 
estimated from the standard deviations of the measurements for 
non-variable stars, are about ay(HB) = 0.13 and ctb(HB) = 0.11 
on the individual frames. 

3. RR LYRAE VARIABLE STARS IN LEO I 

The time-series photometric measurements in the CCD cat- 
alogs were analyzed to detect variable objects and derive their 
light curves. Here we report on the first successful results for 
two CCD's (named #6 and #7 in the ESO WFI Manual) com- 
prising the body of the galaxy (Leo I was centered on CCD #7). 
All results for the remaining CCD's, together with photometric 
measurements and light curves, will be presented in a separate 
paper (G. Clementini et al. 2001, in preparation). Identifica- 
tion of the variables was performed independently on the instru- 
mental v and b data sets, running the program VARFIND (devel- 
oped by Dr. P. Montegriffo) on the photometric catalogs (see 
Clementini et al. 2001 for details). At this stage, the search was 
optimized to find candidate variables near the HB. Stars whose 
standard deviations of the v and b measurements are larger than 
3cr, where a is the rms of non-variable stars at the HB level, 
were flagged as candidate variables. As a result, 96 candidates 
were found, some of which are in the innermost area of Leo I. 
RR Lyrae's allow us to probe the old population all the way to 
the center of the galaxy. 

Period search and definition of the light curve of the sus- 

3 IRAF is distributed by the National Optical Astronomy Observatories, which are 
cooperative agreement with the National Science Foundation. 



pected variables were done using differential photometry in the 
instrumental magnitude system, with respect to some carefully 
chosen, photometrically stable reference star (in general, red 
giants with B-V > 0.6). Period search was performed us- 
ing GRATIS (GRaphical Analyzer of Time Series), a software 
developed at the Bologna Observatory (see Clementini et al. 
2000, 2001). 
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FIG. 1 . — Differential v and b light curves of ab- and c-type RR Lyrae's in 
Leo I from data in CCD #6 (left panel) and #7 (right panel). Different symbols 
are used for data taken in different nights. 

Figure 1 presents the first light curves of RR Lyrae stars in 
Leo I. We have plotted the differential v and b light variations 
for a selection of RR Lyrae stars that sample quite well the typ- 
ical periods. Out of 96 candidate variables, the majority were 
found to be ab-type RR Lyrae's (63 stars), with a few exam- 
ples of c-type RR Lyrae's (11 objects). We also identified one 
probable anomalous Cepheid and 6 candidate short period bi- 
naries, while the remaining 15 objects could not be classified 
unambiguously. Some low amplitude variables, in particular 
c-type RR Lyrae's, may actually have escaped detection. A 
deeper search is in progress using the Alard's (2000) Optimal 
Image Subtraction method (G. Clementini et al. 2001, in prep.). 
Periods, amplitudes, and epochs of maximum light were de- 
rived for a subset of 54 RR Lyrae's (47 RRab's and 7 RRc's) 
with complete light curves, using the GRATIS x 2 Fourier fit- 
ting routine (one harmonic was used for the c-type RR Lyrae's 
and two to five harmonics for the RRab's). The accuracy of 
the derived periods depends on the sampling of the light curves 
and is in general of the order 0.001 d. The amplitude uncer- 
tainties, estimated from the rms of the fit residuals, are ^0.10 
mag. Intensity-average differential magnitudes were obtained 
by integration over the entire pulsation cycle using the best fit- 
ting model curves. By adding the instrumental magnitudes of 
the reference stars, we obtained the mean b, v magnitudes of 
the RR Lyrae's, from which the mean B, V magnitudes were 
calculated via the calibration equations. 

operated by the Association of Universities for Research in Astronomy, Inc., under 



3 



20 

V 

22 



24 



26 



1 1 1 1 1 

_ LEO 1 (CCD n.6) 


1 1 1 | 1 1 1 1 | l 1 


_ «RR Lyrae 




■Cepheid 




1 









1 2 

B-V 



FIG . 2 . — HR diagram of Leo I from a master ALLFRAME combination of all 
data for CCD #6. This CMD is based on the mean calibrated B, V magnitudes 
of stars detected in both passbands. The RR Lyrae variables (filled circles) and 
a Cepheid (filled square) identified in this chip are plotted according to their 
intensity-averaged magnitudes and colors. Lines show the blue and red edges 
of the instability strip of the globular cluster M 3, from Corwin & Carney's 
(2001) photometry, reddened to the value appropriate for Leo I. The red edge 
has been corrected to account for the difference in metallicity with respect to 
Leo I (Walker 1998). 

The location of the detected variables in the HR diagram of 
the outer region of Leo I is shown in Figure 2. The vast majority 
of the data points for the RR Lyrae's are found on the HB in cor- 
respondence to the instability strip, and their mean magnitude 
defines the mean V apparent luminosity of the HB of Leo I. 



4. PULSATIONAL PROPERTIES: CLUES TO THE OLD 
STELLAR POPULATION IN LEO I 

Figure 3 shows the relation between the period and the V 
amplitude for the RR Lyrae's in the present sample. Open 
symbols represent the c-type pulsators (RRc), filled symbols 
are RRab stars. The average period of the Leo I RRab vari- 
ables is (P fl /,)=0?602 (cr=0?059), and the minimum period is 
0?539. The solid lines in Figure 3 represent the Ay-logP 
relations defined by the ab-type variables with clean light 
curves in the globular clusters M3 and M 15 ([Fe/H]=-1.66 
and [Fe/H]=-2.15, respectively, on the Zinn & West's 1984 
scale). These clusters can be considered the prototypes of Oost- 
erhoff type I (Ool) and II (OoII) globular clusters (GCs) in the 
Milky Way (Sandage, Katem, & Sandage 1981; Sandage 1993). 
Also shown are the Ay - log P relations followed by the ab vari- 
ables in M2, w Cen, and NGC6441 ([Fe/H]=-1.62, -1.60, and 
-0.5 dex). The data are from Clement & Shelton (1999) for 
M3, Bingham et al. (1984) and Silbermann & Smith (1995) 
for M15, Lee & Carney (1999a) for M2, Clement & Shelton 
(1999) for oj Cen, and Pritzl et al. (2000) for NGC6441. 
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FIG. 3. — Period-V amplitude relation for the RR Lyrae's in Leo I. Periods 
are in days. Filled and open symbols are RRab and RRc variables, respectively. 
The regression lines show the Ay — log P relations for the ab-type variables 
in the globular clusters M3 and M 15 (solid lines), M2 (dotted line), w Cen 
(dashed line), and NGC 6441 (long dashed line). Dots represent the amplitude- 
period distribution of RR Lyrae variables in Leo II (Siegel & Majewski 2000). 

Figure 3 also shows a comparison with the Ay- log P dis- 
tribution of the RR Lyrae variables in the dSph galaxy Leo II 
(Siegel & Majewski 2000). The pulsational properties of the 
RR Lyrae stars in the two dwarf spheroidals appear quite sim- 
ilar. The characteristics of the ab-type RR Lyrae's qualify 
Leo I as a system intermediate between the Ool and OoII clus- 
ters, similar in this respect to other dwarf galaxies in the Local 
Group (Siegel & Majewski 2000; Cseresnjes 2001; and refer- 
ences therein) and to the LMC (Bono, Caputo, & Stellingw- 
erf 1994). However, both the rather long minimum period of 
RRab's and the presence of RRab's with long periods - large 
amplitudes, indicate a larger similarity with the OoII GCs. 

The implications of such Oosterhoff intermediate properties 
for the old stellar populations in Leo I (and in other dSph's as 
well) rests on the usual interpretation of metallicity as the dom- 
inant parameter in determining the Oosterhoff dichotomy. In- 
deed, OoII clusters are generally very metal-poor while Ool 
clusters are of intermediate metallicity (e.g., Smith 1995). If 
this interpretation is correct, then the distribution of periods and 
amplitudes of Leo I RR Lyrae's implies a metallicity distribu- 
tion for the old population extending from values more metal- 
poor than [Fe/H]=-2.15 dex (the metallicity of the OoII cluster 
M 15), to as metal-rich as [Fe/H]^ -1.6 (the metallicity of the 
Ool cluster M 3). 

Alternative interpretations may exist, though. There are clus- 
ters like M2, w Cen, NGC 6441, and NGC 6388 whose po- 
sitions in the period-amplitude diagram imply that their RR 
Lyrae's are brighter than expected for their metallicity (see Fig- 
ure 3). This observational evidence is interpreted as the Oost- 
erhoff dichotomy being caused by evolution off the zero-age 
horizontal branch (ZAHB), with the Ool variables being on the 
ZAHB and the OoII variables being more evolved (Clement & 
Shelton 1999, Clement & Rowe 2000, Lee & Carney 1999b). 
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Then, a difference in age would exist between the two Oost- 
erhoff groups. Whether this can be a general interpretation of 
the Oosterhoff phenomenon or these are just peculiar clusters 
is still a matter of debate. In this regard, it is interesting to 
note the presence in Leo I of a variable with very long period 
(star 201064 in Figure 3, P = d 889). Although its large ampli- 
tude may be in doubt since the variable falls near the center of 
Leo I where crowding is severe, and is possibly affected from 
straylight from Regulus, variables with long periods have been 
found in NGC 6441 and NGC 6388 (Pritzl et al. 2000). 

A good indication on the mean metallicity of RR Lyrae's 
in Leo I can be derived using the relation derived by Sandage 
(1993) between the average period of the RRab variables 
in the Milky Way GCs and their metallicity, log(P„/ 7 ) = 
-0.092[Fe/H]- 0.389. By applying this relation to the mean pe- 
riod of RRab variables in Leo I, we obtain [Fe/H]^) = -1.82 
for the average metallicity of the old population. Previous val- 
ues of the metallicity of Leo I are quite uncertain, ranging from 
[Fe/H]=-2 to -1 (see Lee et al. 1993 and references therein). 
A recent comparison of the (V — I) color of the red giant branch 
(RGB) of Leo I with the fiducial sequences of Galactic GCs 
yields a mean metallicity [Fe/H]~ -1.8 (Y. Momany et al. 
2001, in preparation). Since the bulk of the red giant stars in 
Leo I have intermediate age (between 2 and 7 Gyr: Caputo et 
al. 1998; Gallart et al. 1999b), this value should be corrected 
(made more metal-rich) by ~ 0.2 dex (cf. Held, Saviane, & Mo- 
many 1999). The fact that the RR Lyrae's are only slightly more 
metal-poor than RGB stars seems to imply that metal enrich- 
ment (from the primordial gas composition to [Fe/H]~ -1.8) 
occurred early in the life of the galaxy, and modest enrichment 
took place in between the early burst that gave rise to the oldest 
population, and the main star formation episode occurred 4 ± 3 
Gyr ago. 

5. THE DISTANCE TO LEO I 

The mean magnitude of the RR Lyrae's in Leo I provides an 
independent method to estimate the distance to this galaxy with 
some degree of confidence. The average apparent luminosity of 
the RR Lyrae's with full coverage of the B and V light curves is 
(V (RR)) = 22.60 ±0.12 mag (standard deviation of the mean, 
48 stars). This dispersion is fully accounted for by the com- 



bined effects of photometric errors, metallicity distribution in 
the RR Lyrae's sample, and evolution off the ZAHB. Given the 
large distance of Leo I, the spread caused by the intrinsic depth 
of the galaxy is insignificant. 

The absolute magnitude of the RR Lyrae's is My (RR) ~ 
0.50 at [Fe/H]= -1.5 in the "long" distance scale, consis- 
tent with (m-M)^ R = 18.53 for the LMC (see Walker 1999, 
and references therein), and 0.75 mag in the "short" scale. 
To correct the absolute magnitude of the RR Lyrae's to 
the metallicity of the Leo I variables, we adopted the RR 
Lyrae luminosity-metallicity dependence given by the rela- 
tion AMy(RR)/A[Fe/H] = 0.2 mag/dex, which is supported 
by both the Baade-Wesselink analysis of Fernley et al. (1998) 
and recent results on the M31 globular clusters (Corsi et al. 
2000). For the metallicity of Leo I we assumed the value de- 
rived in this paper from the mean period of ab-type RR Lyrae's, 
[Fe/H]=-1.82. A reddening E(B-V) = 0.04 ±0.02 mag was 
adopted following Schlegel, Finkbeiner, & Davis (1998). These 
assumptions yield a distance modulus (m-M)o R = 22.04 ±0.14 
mag (<f=256 kpc) in the long scale (and 21.79 mag or 228 kpc 
in the short scale). The quoted error includes the errors on the 
mean RR Lyrae magnitude, the calibration zero point error, and 
the adopted reddening uncertainty. This result is in good agree- 
ment with our recent distance determination from the RGB tip 
method on a different data set (Y. Momany et al. 2001, in 
prep.), yielding (m-M)Jj p = 21.93 ±0.16. The RGB tip method 
is based on the RR Lyrae luminosity-metallicity relation of 
Lee, Demarque, & Zinn (1990), yielding My (RR) = 0.565 at 
[Fe/H]= -1 .5. Our new result is intermediate between previous 
values in the literature (see Lee et al. 1993). 

We thank R Stetson for providing us his set of photometric 
programs, including ALLFRAME. We are grateful to M. Cate- 
lan for helpful discussions regarding the nature of the Oost- 
erhoff dichotomy, to C. Cacciari for useful comments on the 
manuscript, and to an anonymous referee for helpful remarks. 
This study was partially supported by the National projects 
"Processing of large-format astronomical data" (COFIN-98) 
and "Stellar Observables of Cosmological relevance" (COFIN- 
00). 
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